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Abstract. We provide evidence that it is computationally difficult to approxi- 
mate the partition function of the ferromagnetic g-state Potts model when q > 2. 
Specifically we show that the partition function is hard for the complexity class 
^RHLTi under approximation-preserving reducibility. Thus, it is as hard to 
approximate the partition function as it is to find approximate solutions to a 
wide range of counting problems, including that of determining the number of 
independent sets in a bipartite graph. Our proof exploits the first order phase 
transition of the "random cluster" model, which is a probability distribution on 
graphs that is closely related to the g-state Potts model. 



I. Introduction 

Let q be a positive integer. The g-state Potts partition function of a graph G = 
(V, E), with uniform interactions of strength 7 > — 1 along the edges, is defined as 

(1) -Zpotts 

a:V-*[q] e={u,v}&E 

where [q] — {1, . . . , q} is a set of q spins or colours, and 5(s, s') is 1 if s = s', and 
otherwise. The partition function is a sum over "configurations" a which assign 
spins to vertices in all possible ways. Later, we shall widen this definition to allow 
the interaction strength 7 to be a function of the edge e, and allow G to be a 
hypergraph (a generalisation of graph in which edges contain an arbitrary number 
of vertices), but the above restricted case is sufficient for this overview. Mostly 
we shall concentrate in this paper on the ferromagnetic situation, characterised by 
7 > 0. In the ferromagnetic Potts model, configurations o with many adjacent 
like spins make a greater contribution to the partition function Zp otts (G; q, 7) than 
those with few. 

The statistical mechanical model just described was introduced by Potts [25] 
and generalises the classical Ising model from two to q spins. 

Definition ([I]) applies only when q is a positive integer. However, it transpires 
that, regarding q as an indeterminate, ([I]) defines a polynomial in q, and in this 
way we can make sense of the Potts partition function for non-integer q, even 
though the underlying physical model has no meaning. An equivalent, but more 
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concrete way of approaching the partition function when q is non-integer is via the 
Tutte polynomial, which in its "random cluster" formulation is defined as follows: 



where k(V,F) denotes the number of connected components in the graph (V,F). 
The notation is as before, except that now q is an arbitrary real number. Again, 
for simplicity, we are assuming that G is a usual graph (not a hypergraph) and that 
the edge weight 7 is uniform over edges. The multivariate definition can be guessed 
at and will in any case appear later in the paper. For readers who are familiar with 
the classical (x, y)-parameterisation of the Tutte polynomial, the transformation 
between that and the one here is given by 7 = y — 1 and q = (x — l)(y — 1). 

Although ([1]) and (J2J) are formally very different, they define the same polynomial 
in q: see Observation |2j We continue the discussion now in terms of the Tutte 
polynomial 02]), remembering all along that we include as a special case the Potts 
partition function, and as an even more special case that of the Ising model. 
We denote by Tutte^, 7) the computational task of computing Zq^tteiG; q, 7) 
given a graph G as problem instance. Then each pair (g, 7) defines a separate 
computational problem, and we can study the computational complexity of this 
problem as q and 7 vary. It is important to note that q and 7 do not form part 
of the problem instance, which consists simply of the graph G. For the purposes 
of this discussion, we may assume that q and 7 are rational, in order to avoid 
representation issues, but in the main body of the paper we work in the wider 
class of "efficiently approximable" real numbers. 

In a seminal paper, Jaeger, Vertigan and Welsh [18] examined the problem of 
computing Z Tutte (G; q, 7) exactly. In the exact setting, they completely classified 
the complexity of TuTTE(g, 7) for all q, 7 (in fact for all complex q, 7). It transpires 
that TuTTE(g, 7) is #P-hard (i.e., as hard as determining the number of satisfying 
assignments to a CNF Boolean formula), except when q = 1, or when (5,7) is one 
of a finite number of "special points"; in these cases TuTTE(g, 7) is polynomial- 
time computable. 

In light of Jaeger et al.'s strong negative result, attention turned to the question 
of whether Z^^iG] q, 7) could be approximated with arbitrarily small specified 
relative error. In the context of computing partition functions, the appropriate 
notion of efficient approximate computation is the "Fully polynomial randomised 
approximation scheme" or FPRAS, which is rigorously defined in £j5j An early 
positive result was provided by Jerrum and Sinclair [20] , who presented an FPRAS 
for the case q = 2 and 7 > 0, that is to say, for the ferromagnetic Ising model. 
Sadly, no further generally applicable positive results have appeared since then, 
though FPRAS's have been proposed for restricted classes of graphs, e.g., dense 
or degree-bounded [T]. 



(2) 




APPROXIMATING THE POTTS PARTITION FUNCTION 



3 



Greater progress has been made in the negative direction. Goldberg and Jer- 
rum [TJJ showed, under the reasonable complexity-theoretic assumption RP 7^ NP, 
that no FPRAS exists for TuTTE(g, 7) for a wide range of values for the parame- 
ters (5,7). Stated informally, RP 7^ NP is the assumption that there are problems 
in NP that cannot be decided by a polynomial-time randomised algorithm. Intu- 
itively, it is only slightly stronger than the more usual P 7^ NP assumption. As an 
indicative example of what is known, the intractability result of [12] covers the en- 
tire half-plane 7 < —2 except for the tractable case q = 1 and the case q = 2 where 
the problem is equivalent to approximately counting perfect matchings. Similar 
results apply when g/7 < —2. The restriction to planar graphs was treated in a 
follow-up paper [13J. However none of the existing intractability results apply to 
the region q > and 7 > that concerns us here, and which is perhaps the one of 
greatest physical interest. 

Our goal here is to present the first evidence that TuTTE(g, 7) is computa- 
tionally hard in the region q > 2 and 7 > 0, i.e., the region corresponding to the 
ferromagnetic Potts model with q > 2 states. We achieve this, but under a stronger 
complexity-theoretic assumption than RP 7^ NP. To explain this assumption, a 
digression into computational complexity is required. 

The complexity class ^KRUi of counting problems was introduced by Dyer, 
Goldberg, Greenhill and Jerrum [S] as a means to classify a wide class of ap- 
proximate counting problems that were previously of indeterminate computational 
complexity. The problems in #RHQi are those that can be expressed in terms 
of counting the number of models of a logical formula from a certain syntacti- 
cally restricted class. (Although the authors were not aware of it at the time, this 
syntactically restricted class had already been studied under the title "restricted 
Krom SNP" [7]. Yet another terminological variation is to say that problems in 
T^RHIl! enumerate solutions to a linear Datalog program.) The complexity class 
T^RHIl! has a completeness class (with respect to approximation-preserving "AP- 
reductions" ) which includes a wide and ever-increasing range of natural counting 
problems, including: independent sets in a bipartite graph, downsets in a par- 
tial order, configurations in the Widom-Rowlinson model (all [8]), the partition 
function of the ferromagnetic Ising model with mixed external field (i.e., not con- 
sistently favouring one or other spin) [11], and stable matchings [6]. Either all of 
these problems admit an FPRAS (i.e., are efficiently approximable), or none do. 
No FPRAS is known for any of them at the time of writing, despite much effort 
having been expended on finding one. 

All the problems in the completeness class mentioned above are inter-reducible 
via AP-reductions, so any of them could be said to exemplify the completeness 
class. However, mainly for historical reasons, the particular problem f^BIS, of 
counting independent sets in a bipartite graph, tends to be taken as the exemplar 
of the class, much in the same way that Sat has a privileged status in the theory 
on NP-completeness. Our main result is 
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Theorem 1. Suppose that q > 2 and 7 > are efficiently approximable. Then 
#BIS < AP TuTTE(g, 7 ). 

Here, <ap is the symbol for "is AP-reducible to" , and "efficiently approximable" 
is a concept defined in £j5j suffice it to say for now that the rational numbers are 
trivially efficiently approximable. 

One limitation of our result is that the resulting inapproximability of the prob- 
lem TuTTE(g, 7) is conditional on there being no FPRAS for #BIS (and the 
rest of the completeness class), rather than on the weaker assumption NP 7^ RP. 
In fact, we conjecture that t^BIS does not admit an FPRAS. The basis for our 
conjecture is empirical — namely that the collection of known #BIS-equivalent 
problems is growing and that the problem itself has survived its first decade de- 
spite considerable efforts to find an FPRAS. For example, Ge and Stefankovic [9] 
recently proposed an interesting new MCMC algorithm for sampling indepedent 
sets in bipartite graphs. Unfortunately, however, the relevant Markov chain mixes 
slowly [10] so even this interesting new idea does not give an FPRAS. 

Despite the fact that our results are limited by a strong complexity-theoretic 
assumption, we feel there are counterbalancing strengths that fully justify this 
investigation. One is the range and intrinsic interest of the problem under consid- 
eration. Whether in the guise of the Potts partition function, or of the Tutte plane, 
the computational complexity of TuTTE(g, 7) has received considerable attention 
since it was first studied by Jaeger et al. [IS] : see, for example, [JJ EH [28j [29]. So 
it seems worth striving for a complexity classification even under a strong assump- 
tion such as the one we are making. The situation is similar to working with the 
Unique Games Conjecture in the area of approximation algorithms for optimisa- 
tion problems, or employing the class PPAD in analysing the complexity of Nash 
equilibria. Futhermore, Theorem [TJ has a wide range of applicability, covering as it 
does the whole region q > 2, 7 > 0, which, in the classical parameterisation of the 
Tutte polynomial, equates to the entire upper quadrant of the Tutte plane above 
the positive branch of the hyperbola if 2 = {(%, y) '■ — l){y — 1) = 2}. Note that 
the #BIS-hard region extends right to the tractable hyperbola Hi- 

Another potential strength of the work is that the reduction introduces a novel 
technique that may have wider applicability. The idea is conceptually simple 
and can be sketched informally here. In the first step, we reduce ^BIS to a 
hypergraph version of the Tutte polynomial. (The conventional Tutte polynomial 
may be recovered as the specialisation to 2-uniform hypergraphs.) This step, if not 
routine, is at least standard in its techniques. After this, we show how to simulate 
each hyperedge containing t vertices by a graph gadget with t distinguished vertices 
or terminals. 

At this point we exploit the first order phase transition that is a feature of the so- 
called random cluster model when q > 2. The configurations of the random cluster 
model on a graph G are spanning subgraphs of G, which are weighted according 
to the numbers of edges and connected components they contain. As formulated 
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in (j2J), the Tutte polynomial is the partition function of this model. The gadget is 
designed and carefully tuned so that is has two coexisting "phases" : one in which 
the random cluster configurations (spanning subgraphs) have a large connected 
(or "giant") component, and one in which they don't. We show that it is possible 
to arrange for the t terminals to be (with high probability) in a single component 
in one phase and in t distinct components in the other. This provides us with a 
bistable gadget that simulates a potentially large hyperedge using many 2-vertex 
edges. Note that AP-reductions often exploit phase transitions, playing one class 
of configurations off against another. See the examples in [5] and [23]. What is 
new here is that, as far as we are aware, this is the first time anyone managed to 
derive stronger complexity results using the complex phase transitions that arise in 
actual models studied in statistical physics. Unfortunately, the delicate nature of 
the gadgets needed to exploit this kind of phase transition does lead to significant 
technical complexity in our analysis. 

Non-trivial phase transitions have been used in the past to rule out certain 
natural Markov chain Monte Carlo (MCMC) approaches to approximate counting 
problems. Indeed, the first-order phase transition in the random cluster model 
was already exploited by Gore and Jerrum [14J to demonstrate that the Swendsen- 
Wang algorithm is not always effective for the ferromagnetic Potts model. (See also 
Borgs et al. [I] for a more thorough working out of this idea.) In those applications 
the aim was to rule out a certain algorithmic approach, namely MCMC, whereas 
here our goal is to show inherent intractability. 

Finally, note that Theorem [1] establishes #BIS-hardness of TuTTE(g, 7) but not 
#BIS-equivalence. It would be very interesting to know whether there is an AP- 
reduction from TuTTE(g, 7) to ^BIS. Note that the complexity of approximate 
counting is complicated. Bordewich [3] has shown that if any problem in #P fails 
to have an FPRAS, then there is an infinite approximation hierarchy within #P. 

The rest of the paper is laid out as follows. Sections [2] and [3J introduce the 
models we shall be working with. Section H] is the technical heart of the paper. 
It analyses the random cluster model on some families of graphs, building up 
to the constuction of the hyperedge simulation gadget. Sections [5] and [6] cover 
the computational framework we work in, while ^71-fTul present and analyse the 
reductions themselves. The final £ fTT1 turns the hypergraph intractability result on 
its head by focusing on the 3-uniform case. By reduction to the graphic case, we 
obtain an FPRAS for the partition function. 

2. The Tutte polynomial of a hypergraph and the Potts partition 

function 

Let H = (V, E) be a hypergraph with vertex set V and hyperedge (multi)set £. 
Following the usual convention for the Tutte polynomial [26], a hypergraph is 
allowed to have parallel edges. This is why £ is a multiset. Also, the vertices in a 
particular hyperedge need not be distinct (so each hyperedge can be viewed as a 



G 



LESLIE ANN GOLDBERG AND MARK JERRUM 



multiset). The multivariate Tutte polynomial of H is defined as follows 

Z Tutte (i/;g,7) = ^g K{V - F) n^ 

Fee feT 

where q and 7 = {if} fee are commuting indeterminates and k(V, J 7 ) denotes the 
number of connected components in the subhypergraph (V, J 7 ). (Two vertices u, v, 
are in the same component of (V, J 7 ) if u = v, or there is a sequence /i,...,/<G J 7 
of hyperedges with u 6 fx, v e fi and fi fl 7^ for 1 < i < £.) This 
partition function was studied (under a different name) by Grimmett [15J. An 
undirected graph G can be viewed as a 2-uniform hypergraph (a hypergraph in 
which every hyperedge has size 2). In this case, Z^atteiG; q,~f) coincides with the 
usual definition of the multivariate Tutte polynomial [26J. 

Let q be a positive integer. The (/-state Potts partition function of H is defined 
as follows: 

Z Potts (H;q n )= J2 II C 1 + I «€/})), 

<r:V^[q] }&£ 

where [q] — {1, . . . , q} is a set of q spins or colours, and S(S) is 1 if its argument 
is a singleton and otherwise. The partition function is a sum ranging over 
objects, in this case assignments of spins to vertices, which are often referred to 
as "configurations". 

The following observation is due to Fortuin and Kastelyn. 

Observation 2. If q is a positive integer then Z Potts (H] q, 7) = Z^^H; q, 7). 

Proof. A proof can be found in [51 Proposition 3.1] or in [15] . but we include it 
here for completeness. The argument is a straightforward generalisation of the 
standard proof for the graph case. See [26, Theorem 2.3]. Consider 

E EIlMW^he/)). 

Now "integrate out" (i.e., explicitly sum over) T to see that this quantity is 
equal to Z PottB (H; q, 7) and integrate out o to see that this quantity is equal to 
Z Tuttc (H;q,j). □ 

The ferromagnetic Potts model corresponds to the Potts model in the special 
case in which the edge weights 7/ are non-negative. In this case, a mono-chromatic 
edge contributes more weight than an edge with multiple spins. 

For a subset J 7 C £ of the hyperedges of a graph, we use 7(J r ) to denote YlfeT^f- 
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3. The Random Cluster Model 

Consider a graph G = (V,E). Every edge e G E is associated with a quantity 
p(e) G [0, 1]. Then for a set of edges AGE define 

P(G;A,q,p) = q^l[p(e) J] (l-p(e)). 

Let 

Z rc (G; g,p) = ]T P(C7; A, q,p) = Z Tutte (G; q, 7 ) I^ 1 ~ p(e)), 

AC_B ee£ 

where 7 e = p(e)/(l — p(e)). Then the probability of edge-set A in the random 
cluster model is given by 

Z IC (G;q,p) 

The random cluster model refers to the distribution KC(G;q,p), in which a 
subset A of edges is chosen with probability P(G;A,q,p). The difference between 
^Tutte an d Z rc is simply one of parameterisation. Nevertheless, the change of 
parameter is useful, as it allows us to employ probabilistic terminology and exploit 
existing results from the random graph literature. 

A central ingredient which helps us to understand the random cluster model is to 
compare it to the (multivariate) Erdos-Renyi model of a random graph. Consider 
a graph G = (V, E). Every edge e G E is associated with a quantity p'(e) G [0, 1]. 
In the Erdos-Renyi model ER(G; p'), a subset A of edges is chosen with probability 

Hp'(e) n (l-p'(e)). 

e£A eeE\A 

Thus, to choose a configuration A in the ER(G;p') model, each edge e G E is 
included in A independently with probability p'(e). 

3.1. Stochastic domination results. We will use several simple stochastic dom- 
ination results, which are close in spirit, and indeed in their proofs, to results of 
Holley [IB] : see, in particular, Theorem (6) of that article and its proof. 

If A + and A~ are disjoint subsets of E, let RC(G; q,p; A + , A~) be the random 
cluster model conditioned on the fact that the chosen subset A contains every edge 
in A + and no edges in A~ . (To avoid trivialities, we assume that no edge e G A + has 
p(e) = and that no edge e G A~ has p(e) = 1.) Similarly, let ER(G;p'; A + , A~) 
be the Erdos-Renyi model with this conditioning. We give conditioned versions of 
the stochastic domination results. 

Lemma 3. Consider a graph G = (V, E) in which each edge e G E is associated 
with a quantity p(e) G [0,1]. Let A + and A" be disjoint subsets of E such that 
p(e) > for e G A + and p(e) < 1 for e G A~ . Suppose q > 1. Then the random 
cluster model RC(G; q,p; A + , A~) is stochastically dominated by the Erdos-Renyi 
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model ER(G;p; A + , A ) in the sense that we can select a pair (A, A') such that 
A is drawn from RC(G; q,p; A + , A~) and A' is drawn from ER(G;p; A + , A') and 
AC A'. 

Proof. Let Ai be the "heat-bath on an edge" Markov chain for RC(G; q,p; A + , A~) 
and let Ai' be the "heat-bath on an edge" Markov chain for ER(G; p\ A + , A~). 
Start Ai in state A and Ai' in state A' Q where A C A' and A and A' Q contain 
every edge in A + and no edges in A~ . We can then couple the evolution of the 
chains from the i'th pair of states (Ai, A[) as follows, guaranteeing that A i+ i C 
A' i+1 . Choose the same edge e G E\ A + U A~ in both chains. In Ai the probability 
of putting e in A i+ i is either p(e) or p(e)/(p(e) + q(l — p(e))), each of which is at 
most p(e), which is the probability of putting it in A' i+1 . □ 

Lemma 4. Consider a graph G = (V, E) in which each edge e G E is associated 
with a quantity p(e) G [0, 1]. For each edge e E E, let p'(e) = p(e)/q. Let A + and 
A~ be disjoint subsets of E such thatp(e) > fore G A + and p(e) < 1 fore G A~ . 
Suppose q > 1. Then the Erdos-Renyi model ER(G; p'\ A + , A~) is stochastically 
dominated by the random cluster model RC(G; q,p; A + , A~) in the sense that we 
can select a pair (A, A') such that A is drawn from ER(G; p'; A + , A~) and A' is 
drawn from RC(G; q,p; A + , A~) and A C A' . 

Proof. Let M. be the "heat-bath on an edge" Markov chain for ER(G;p'; A + , A~) 
and let Ai' be the "heat-bath on an edge" Markov chain for RC(G; q,p; A + , A~) 
Start M in state A and M' in state A' where A C A' and A and A' contain 
every edge in A + and no edges in A~ . We can then couple the evolution of the 
chains from the i'th pair of states (Ai, A^) as follows, guaranteeing that A i+ i C 
A' i+1 . Choose the same edge e G E\A + UA~ in both chains. In Ai' the probability 
of putting e in A' i+1 is either p(e) or p(e)/(p(e) + q{l — p(e))), each of which is at 
least p'(e), which is the probability of putting it in A i+ ±. □ 

Lemma 5. Consider a graph G = (V, E) in which each edge e G E is associated 
with a quantity p(e) G [0, 1]. Suppose that, for each edge e G E, p'(e) > p(e). Let 
A + and A~ be disjoint subsets of E such that p(e) > for e G A + and p(e) < 1 
for e G A~ . Suppose q > 1. Then the random cluster model RC(G; q,p; A + , A~) is 
stochastically dominated by the random cluster model RC(G; q,p'; A + , A~) in the 
sense that we can select a pair (A, A') such that A is drawn from the distribution 
RC(G; q,p; A + , A~) and A' is drawn from RC(G; q,p'; A + , A') and A C A'. 

Proof. Let Ai be the "heat-bath on an edge" Markov chain for RC(G; q, p; A + , A~) 
and let Ai' be the "heat-bath on an edge" Markov chain for RC(G; q, p'; A + , A~) 
Start Ai in state Aq and Ai' in state A' Q where Aq C A' and Aq and A' Q contain 
every edge in A + and no edges in A~ . We can then couple the evolution of the 
chains from the i'th pair of states (Ai, A-) as follows, guaranteeing that A i+ i C 
A' i+1 . Choose the same edge e G E\A + UA~ in both chains. If the endpoints of the 
edge are in the same component of Ai — e then they are in the same component of 
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A' { — e. The probability of putting e in A i+1 is p(e) and the probability of putting 
e in A' i+l is p'(e), which is at least as big. Otherwise, the probability of putting e 
in A i+1 is p(e)/(p(e) + q(l — p(e))). This is at most p'(e)/(p'(e) + g(l — p'(e))) and 
the probability of putting e in A' i+1 is at least p'(e)/(p'(e) + q(l — p\e))). □ 

3.2. The Fundamental Lemma of Bollobas, Grimmett and Jansen. An- 
other way to compare the random cluster model to the Erdos-Renyi model is to 
use the following lemma, which is a multivariate version of Bollobas, Grimmett 
and Jansen's "Fundamental Lemma", [2J Lemma 3.1]. To do the comparison, the 
lemma should be applied with r = 1/q so that the distribution RC(G[Vi]; rq,p) in 
the statement of the lemma below is identical to the distribution ER(C7[Vi]; p). 

Let < r < 1 be fixed. Given a subset A of edges chosen from RC(G; q,p), 
colour each component of (V, A) either red, with probability r, or green, with 
probability 1 — r: different components are coloured independently of one another. 
The union of the red components is the red subgraph and the union of the green 
components is the green subgraph. Let R be the set of vertices in the red subgraph. 
Given a subset U of vertices, let G[U] be the subgraph of G induced by U. 

Lemma 6. (Bollobas, Grimmett and Jansen [21 Lemma 3.1].) Let V\ C V. Con- 
ditioned on R = V\, the red subgraph is distributed according to KC(G[Vi];rq,p) 
and the green subgraph is distributed according to RC(G[V \ Vi]; (1 — r)q,p). Con- 
ditioned on R — V\, the red subgraph and green subgraph are independent of each 
other. 

Proof. Let V 2 = V\V 1 . Let E l = E n V} 2) and E 2 = E n V 2 {2) , where the notation 

(2) 

Vi denotes the set of unordered pairs of vertices from V\. Now consider A% C E% 
and A 2 C E 2 . The (unconditional) probability that the red subgraph is (Vi,Ai) 
and the green subgraph is (V 2 , A 2 ) is 

P(G; At U A 2 , q, p) r< v ^\l - r ) K(V2 ' A2) 
Z TC (G;q,p) 

so the conditional probability, conditioned on R = V\ is 

P(G; A\ U A 2 , q, p) r < v ^\\ - r ) Kiy2 ' M) 
E A[ c Ei ,a' 2 ce 2 P(G; A[ U A' 2 , q,p) r"(*.Ai)(i - 

which is 

P{G[V X \; A u rq,p)P(G[V \ VA; A 2 , (1 - r)q,p). 

□ 

4. The random cluster model on some natural graphs 

In this section we consider the random cluster model on a clique, and also on 
a pair of connected cliques. The latter is used as a gadget in our constructions. 
First, however, we need a technical lemma, 
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4.1. A technical lemma. The following lemma is not about the random cluster 
model, but we will use it in our analysis of the random cluster model on a pair of 
connected cliques. 

Lemma 7. Suppose we have a partition of the set [u] into s blocks of size at 
most f max - Randomly colour a subset of the elements of [v\ yellow according to 
a Bernoulli process with success probability £. Independently, randomly colour a 
subset of the elements of [u] blue according to a Bernoulli process with success 
probability £ (so an element can be coloured yellow, blue, both or neither). We 
say that a block of the partition is bicoloured if it contains both yellow and blue 
elements. Then 

Pr(no block is bicoloured) < [(1 - i) u/s {2 - (1 - £) u/s )] s 

(which is increasing with s and decreasing with u) and 

Pr(some block is bicoloured) < v[l - (1 - £)^—] 2 

(which is increasing in both v and u majX ). 

Proof. Let the block sizes be u±, . . . , v s , so that max.,- Vj < u nmx . Observe that 

Pr(j'th block contains a yellow element) = 1 — (1 — ^) Uj = (say), 

and that the same bound applies, of course, to blue elements. Thus, using the 
convention that the index j always ranges over 1 < j ' < s, 

Pr(some block is bicoloured) < < - (1 - 0" max ] 2 , 

3 3 

which gives the second bound. Clearly, this bound is monotonically increasing in 
v and z/ max as claimed. 
Now for the first bound. 

Pr(no block is bicoloured) = - f ?) = - Q + Q 

3 3 3 

=ri( i -^n( i +o) 

3 3 

(3) =(i-o u U( 1+ ^- 

3 

To get an upper bound we are interested in evaluating the supremum of 

(4) n( i +o)=ri[ 2 -( 1 -^] 

3 3 

over the domain defined by the linear inequalities < Vj < v and . Vj = v. We 
consider this as an optimisation problem over W even though the Uj are all integers; 
of course, this will if anything only increase the supremum. We are considering a 
continuous function over a closed, bounded set, so the supremum is achieved at 
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some point; we claim that this (unique) point is V\ — ■ ■ ■ — v s — vjs. For if not, 
then at least one pair, say v\ and z/ 2 would be unequal. But then it is easily checked 
that replacing v x and v 2 by their average would increase [2 — (1— ^)' /1 ][2 — (1— £)" 2 ], 
and hence increase the right hand side of (jl]): a contradiction. Substituting uj = 
vjs into equations (j4j) and then d^J), we obtain 

(5) Pr(no block is bicoloured) < (1 - i) u [2 - (1 - £) v/s ] s 

(6) = [(i - er /a (2 - (i - er /s )r> 

as desired. 

In only remains to verify the monotonicity claims about (JHD- Let u(u) = (1 — 
g (0,1]. Regarding s > as fixed, u{y) decreases monotonically with v. 
Also, [u(2 — u)} s increases monotonically as a function of u in the range (0, 1). 
Thus, expression decreases monotonically with v. 

Now regard v as fixed and make the change of variable s = ax, where a = 
z/ln((l— O 1 ); an d note that x > 0. Then expression §5§ becomes [e _1 (2— e - l / x Y] a . 
Thus, it is enough to show that fix) = xln(2 — e~ l ^ x ) increases monotonically 
with x. Now 

f"(x) = -2x- 3 e- 1/x (2 - e- 1/x y 2 < 0, 

and 

f'{x) = ln(2 - e~ 1/x ) - x- 1 e~ 1/x {2 - e^)" 1 ^0 as x -> oo. 
These two facts imply f'(x) > for x > 0. □ 

4.2. The random cluster model on a clique. Bollobas, Grimmett and Jansen 
[2] studied the random cluster model on the complete N- vertex graph K^. More 
detailed analyses have since been performed, for example by Luczak and Luczak 
|24J, but the approach of the earlier paper is easier to adapt to our needs. 

For fixed q and a fixed constant A, they studied the distribution RC(KN,q,p) 
where p is the constant function which assigns every edge e of the value 
p(e) = X/N. They show that there is a critical value A c , depending on q, so that, 
if A > A c then, as iV — > oo, with high probability a configuration A drawn from 
KC(Kn, q,p) will have a large component (of size linear in N) and otherwise, with 
high probability the largest component will be much smaller (of size logarithmic 
in N). 

For q > 2, the critical value A c is defined as follows. 

It is important for our analysis that A c < q (see [2, p. 16], or by calculus). We 
define S = (q - A c )/2 > and A = A c + 5. Let = (q- 2)/(q - 1). 

We will use the following lemma, which follows from Theorem 2.2 and Equa- 
tion (5) of 0. 
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Lemma 8. Fix q > 2 and define X and 6 as above. Let p be the constant function 
which assigns every edge e of K N the value p(e) = X/N. Let A be drawn from 
KG(K N ,q,p). The probability that A has a connected component of size at least 
9N tends to 1 as N — )■ oo. 

4.3. The random cluster model on a pair of connected cliques. Let r be 

the complete graph with vertex set Vp = K U T. Let Ep denote the edge set 
of r and let N — \K\ and t = \T\. Let denote the set of unordered pairs of 
distinct elements in K and define similarly. Let g be a value in [0, 1]. Define 
p as follows. 



deleted) as a gadget to simulate the contribution of a hyperedge on the set T to the 
multivariate Tutte polynomial of r. Thus, we refer to vertices in T as "terminals" 
of the graph r. For a subset A C Ep, let Y(A) be the number of connected 
components in the graph (Vp, A \ T^) that contain terminals. 

A remark about the gadget r and its eventual use. When we come to use the 
gadget, the edges in will not be present. It is for this reason that we are 
interested in the structure of connected components in r in the absence of these 
edges, and specifically the random variable Y(A). However, it turns out that 
the key properties of the gadget are easier to verify if we work with a random 
cluster distribution associated with r, exactly as given above, with the edges 
present. Informally, the appropriate "boundary condition" for the gadget is the 
one in which the terminals are joined with probability 1. 

The following two lemmas establish some useful properties of the gadget. The 
second shows that distribution of Y(A) is concentrated at the extremes of its 
range, i.e., Y = 1 or Y = t. This concentration property holds for a wide range 
of values for the edge probability g. The first lemma, which is easier, shows that 
we can tune g so that the balance of probability between those extremes is the 
desired value 7. Later, in Lemmas [12] and [13j we shall show that a sufficiently 
close approximation to this g can be efficiently computed. 

Lemma 9. Fix q > 2 and let X = X c + (q — X c )/2. Fix a weight 7 > and let No be 
a sufficiently large quantity depending on q and 7. Suppose a number of terminals 
t > 1 is given and fix N > max{£ 16 , A^ }. Then there is a parameter g satisfying 
N~ 3 < g < X/N < I such that, if A is drawn from RC(i~'; q,p) then 



Proof. Define 9 (depending on q) as in Section 14.21 Let ip(g) = Pt(Y(A) = 
t)/ Pt(Y(A) = 1). We will use stochastic domination to show 




(7) 



Pr(Y(A) = 1) = 7 Pr(y(A) = t). 
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• i/j(q) is monotonically decreasing as a function of g, 

• ij(N- 3 ) > I/7, and 

• ip(X/N) < 1/7. 

Since ip(g) is a rational function in g and the denominator is never zero, we conclude 
that ip(g) is continuous in g and there is a value g G (N~ 3 , X/N) such that ip(g) — 
I/7. This gives ([7]). Note that the lower bound N~ 3 for g is very crude, but this 
is all that we will need. 

First, to show that ip(g) is monotonically decreasing as a function of g, note from 
Lemma [5] that Pr(Y(A) = t) is monotonically decreasing in g and Pr(Y(A) = 1) 
is monotonically increasing in g. 

Next, to show that tp(N~ 3 ) > I/7, we will assume g = N~ 3 (and that iV is 
sufficiently large) and we will show Pr(Y(A) = t) > 1/(1 + 7), which suffices. 
By Lemma [3j Pr(Y(A) = t) > Pt(Y(A) = t) where A is drawn from ER(.T;p). 
Now the probability that A n = is at least 1 - (^) N~ 3 . In this case, the 
probability that a particular pair of terminals is connected in A \ is at most 
iV x iV _3//4 x jV~ 3 / 4 = iV" 1 / 2 , and the probability that there exists a connected pair 
is at most tfN- 1 / 2 < N' 3 '\ So Pr(F(A) = t)> 1 - (f)N~ 3 - N~^ 8 > 1/(1 + 7). 

To finish the proof of (ED, we will show that ip(X/N) < I/7. To do this, we will 
assume g = X/N (and that iV is sufficiently large) and we will show Pt(Y(A) = 
1) > 7/ (I + 7), which suffices. We will again use stochastic domination to compare 
the random cluster model to the Erdos-Renyi model, but this time we need some 
conditioning. 

First, we will show that the probability that the graph (K, A n K^) has a 
connected component of size at least 6W tends to 1 as iV — )■ 00. To do this, let A* 
be drawn from the distribution RC(r,g,p) where 

3( e ) = / *?> if e G K^, 
\ 0, otherwise. 

Since p(e) < p(e), Lemma \5\ guarantees that the probability in question is at least 
the probability that the graph (K, A*) has a connected component of size at least 
#iV. By Lemma El this probability tends to 1 as N — > 00. 

Next, we will consider the generation of configuration A from the distribution 
RC(.T; q, p) as follows. First, we will select a set A + C from the appropriate 
induced distribution. Letting A~ = K^\A + , we will select A from the distribution 
RC(-T; q,p; A + , A~). We will finish by showing that, as long as (K, A + ) has a 
connected component of size at least ON, Pt(Y(A) = 1 | A + C A, A n A~ = 0) is 
greater than 7/(1+7). To do this, let p'(e) = p(e) / q and let A be a random variable 
drawn from ER(r;p'; A+, A~). By Lemma II Pr(Y(A) = 1 | A + C A, A n A~ = 
0) > P T {Y{A) = 1). Now in A, the probability that a particular terminal is not 
connected to the large component is at most (1 — q~ 1 N~ 3 /' i ) 9N < exp(— Oq 
and the probability that there exists a terminal that is not connected to this 
component is at most t exp(— Oq^N 1 ' 4 ). So Pr(Y(A) = 1) > 7/(1 + 7). □ 
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Lemma 10. Fix q > 2 and let A = A c +(g — A c )/2. Fix a weight 7 > and let N be 

a sufficiently large quantity depending on q and 7. Suppose a number of terminals 
t > 1 and a tolerance < 77 < 1 are given and fix N > max{r; 16 , r/ _1//8 , A^ }- For 
every value of g in the range [N~ 3 , \/N], if A is drawn from RC(-T; q,p) then 

(8) Pr(l < Y(A) <t) < 77. 

Proof. Let C±, C2, ... be the connected components of A D K^ 2 \ ordered in non- 
increasing size. We are going to be relying on the phase transition of the random 
cluster model. The main fact that we will use is that \C\\ is likely to either be 
very small (around order log(iV)) or very large (a constant fraction of N). Thus, 
it is unlikely that \C\\ is close to iV 1//s . We will not need much detail about the 
phase transition. We will show 

(9) Pr ((|d| < iV 1/s ) A (Y(A) < t)) < rj/2, 
and 

(10) Pr ((Idl > iV 1/s ) A iY(A) > 1)) < 77/2. 

(Actually, the event mentioned in (flOi) holds with all but exponentially small prob- 
ability.) The required inequality jS]) follows directly from ([9]) and (|T0|) . 

Inequality (Q is easier. To generate a configuration A from RC(.T; q,p) we first 
select a set A + C from the appropriate induced distribution. Then, letting 
A~ = \ A + , we select A from the distribution RC(.T; q,p; A + , A~). The size 
of C\ is entirely determined by A + . If A + has a component of size greater than 
iV 1 / 8 then 

Pr((|d| < A^ 1/8 ) A (Y(A) < t) I A + C A, An A~ = 0) = 0. 

Otherwise, this probability is equal to Pr(Y(A) < t \ A + C A, A (1 A~ = 
0). By Lemma [3j this is at most Pr(Y(A) < t) where A is generated from 
ER(r;p;A+,A-). _ 

So now consider A. Fix attention on two terminals, and colour the vertices in 
K that are adjacent to the first terminal yellow, and those adjacent to the second 
terminal blue. We are in the situation of Lemma [7] with v — N, z/ max < iV 1//8 
and £ = iV~ 3 / 4 . Thus the probability that there exists a bicoloured block (i.e., 
that the terminals are connected via some connected component Cf) is at most 
N[l - (1 - jV-3/4)ivV8j 2 < N p__ 1 + jy-s/sp = JV -i/4_ Thug the probability that 

there exists a pair of connected terminals is at most ( 2 )iV _1 / 4 < N~ 1 / 8 /2 < rj/2. 
Note that the event "there exists a pair of connected terminals" is the same as the 
event Y(A) <t. 
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At a high level, the path to establishing inequality ([TO]) is as follows. We define 
events £ ( "large component" ) and W ( "weighty components" ) and establish 

(11) Pr ((|Ci| > A^ 1/8 ) A < 77/6, 

(12) Pr(£ A nW) < 77/6, and 

(13) Pr (W A (Y(A) > 1)) < 77/6. 

Then the required inequality ffTD]) follows by elementary algebra of sets (events). 

We will say that an event holds "with high probability" (abbreviated whp) if 
the probability that it holds is o{l/Q{n)) for any fixed polynomial Q. 

First, we prove inequality (TTTj) . Let Si,S%,... be the connected components of 
A and let Ci, C2, . . . be the sets Si D K^ 2 \ S2 H K^ 2 \ . . ., ordered in non-increasing 
size. (Thus, {Cj} is a coarsening of {Cj}.) Event £ is |Ci| > iV 15 /i6. Let g 1 be the 
event (|Ci| > N 1 ^ 8 ) A We are interested showing that the event £\ is unlikely. 

Construct the red subgraph and green subgraph of (Vp, A) as in Lemma M with 
r = l/q. Let £2 be the event 

(\R HK\< N/q + iV 15/16 ) A (Ci Q R) A > iV 1/8 ). 

We will show 

Pr(£ 2 ) > Pr(^)/g 2 

and 

Pr(£ 2 ) < r//(6g 2 ), 

which proves (II ip . The first of these follows from Pr(£ 2 | £\) > q~ 2 , which 
follows since Pr(Ci C R \ £1) = q^ 1 (by the definition of the red subgraph) and 
Pt(\R r)K\< N/q + N 15 / 16 I £ x A (Ci C R)) > q- 1 (also by the definition of the 
red subgraph — with probability q~ 1 , red is the rarest colour among vertices in 
K \ Ci where Cj is the element of {Cj} containing C\). 
Now the probability of event £2 is at most 

Pr ((d C R) A (|d| > A^ 1/8 ) I \R (1 K\ < N/q + iV 15 / 16 ) 

which is at most the conditional probability that A[R D K\ has a component of 
size at least iV 1 / 8 , conditioned on |i? D K\ < N/q + N 15 / 16 . But, by Lemma El 
A[i?ni( r ] is an Erdos-Renyi random graph with |i2DiC| < (g _1 + A^ _1 / 16 )A^ vertices 
and subcritical edge probability g < X/N < \R D K^ 1 . So the probability that 
it has a component of size at least N 1 ^ 8 is at most e(N), where e(-) is smaller 
than any inverse polynomial [121 Proof of Theorem 5.4]. Thus, we have shown 
P r (^2) < v/i^Q 2 ) (provided N is sufficiently large), so we have proved (TTTj) . 

Next we prove inequality (]TJ). The event W is that \C\\ H h |C S | > N 15 ^ 16 , 

for some s < 2N 5 ^ 16 . The component C\ is in general composed of a number of 
components from {Cj}. However, the number of constituent components cannot 
be larger than the number \An (K x T) \ of edges joining K to T. This number can 
be upper-bounded using a bounding configuration A drawn from ER(i~';p). We 
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see from a Chernoff bound that whp the number \A R (K x T) of edges joining 
K to T is less than 2tiV 1 / 4 < 2iV 5 / 16 , i.e., twice the expected number, so whp, at 
most 2iV 5 / 16 components from {Cj} contribute to C\. Thus Pr(W | C) = 1 — e(N) 
and (H2D follows. 

Finally we prove inequality ({13]) . As before, generate A by selecting a set A + C 
f^( 2 ) from the appropriate induced distribution, letting A~ = K^ 2 > \ A + , and 
selecting A from the distribution RC(-T; q,p; A + , A~). We will condition on the 
event W (which is entirely determined by A + ). As above, 

Pr (Y(A) = 1 \ A + C A,Af) A~ =$) > Pv(Y(A) = 1), 

where p'(e) = p(e)/q and A is a random variable drawn from ER(i~'; p'\ A + , A~). 
Consider two terminals i,j in T. We can use Lemma [7J to find an upper bound 
for the probability that i and j are not connected in A via one of the components 
Ci, . . . , C s . In A, edges from T to K are selected independently with probability 
g-ijV -3 / 4 . If we colour vertices in K adjacent to i (respectively, j) yellow (respec- 
tively, blue) then we are in the situation of Lemma [TJ with £ = q~ l N^^. From 
the remarks about monotonicity in s and z/, we may assume for an upper bound 
that v = iV 15//16 and s = 2iV 5 / 16 . The probability we want to bound is that of not 
having a bicoloured component. Observe 

(1 - Ws = (1 - q-'N-V 4 )^ 8 < 1 - ig' 1 ^- 1 / 8 , 

for iV sufficiently large, and since the function u i— >■ u(2 — u) is increasing in the 
range (0, 1), 

(1 - 0^(2 - (1 - 0^) < (1 - \q- x N-^){\ + ^N- 1 / 8 ). 

Applying Lemma [71 

Pr(* + j in A) = (1 - ±q- 2 N^ A ) 2N5/16 < eX p(-|g- 2 iV 1 /i6 ) . 

So whp i ~ j. It follows that whp all vertices in T are connected to each other, and 
so Yr(Y(A) — 1) > 1 — s(N). This deals with (fT3|) and completes the proof. □ 

4.4. The random cluster model on a clique connected to an independent 
set. Construct r as in Section I4T51 Let r = (Vp, E r \ T^) be the graph derived 
from r by deleting edges within T. The vertices in T are called the "terminals" 
of r. Let 7 = {7 e } be the set of edge weights defined by 7 e = p(e)/(l — p(e)). 

For an edge subset A' C E r \T^ 2 \ let k'(Vt, A') denote the number of connected 
components that do not contain terminals in the graph (Vp, A'). Let A k denote the 
set of edge subsets A' C Ep\T^ for which the terminals of (Vp, A') are contained 
in exactly k connected components. Let A = IJfcg[i] A k (this is the set of all edge 

subsets of r). Let Z k be q~ k times the contribution to Z^ utte (r; q,^f) from edge 
subsets A' E A k . Formally, Z k = J2 A , eAk q K,{yA 'hW- Let Z = Y! k =i zk - 

We will use the following lemma to apply Lemmas [9] and [10] from Section 14.31 in 
our reductions. 
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Lemma 11. Z k jZ = Pt(Y(A) = k), where A is drawn from RC(-T; q,p). 
Proof. From the definitions of the random cluster model, 



5. Computational problems, fully polynomial randomised 
approximation schemes and efficiently approximable real 

NUMBERS 

Fix real numbers q > 2 and 7 > and consider the following computational 
problem, which is parameterised by q and 7. 
Problem: TUTTE(g, 7). 
Instance: Graph G = (V,E). 

Output: Z Tuttc (G; q, 7), where 7 is the constant function with 7 e = 7 for 
every e G E. 

We will have much more to say about computational approximations of real 
numbers below. For the moment, it may help the reader to think of q and 7 as 
being rational. Jaeger, Vertigan and Welsh [18] have shown that TuTTE(g, 7) is 
#P-hard for every fixed q > 2 and 7 > 0. Thus, it is unlikely that there is a 
polynomial-time algorithm for exactly solving this problem. (If there were such 
an algorithm, this would entail P = $P, and of course P = NP.) 

We are interested in the complexity of approximately solving TuTTE(g, 7). We 
start by defining the relevant concepts. A randomised approximation scheme is an 
algorithm for approximately computing the value of a function / : X* — > M.. The 
approximation scheme has a parameter e > which specifies the error tolerance. 
A randomised approximation scheme for / is a randomised algorithm that takes 
as input an instance x G S* (e.g., for the problem TUTTE(g, 7), the input would 
be an encoding of a graph G) and a rational error tolerance e > 0, and outputs a 
rational number z (a random variable of the "coin tosses" made by the algorithm) 
such that, for every instance x, 





which is what we require, once we cancel a factor of 



q H (!-P(e)) 



from the numerator and denominator. 



□ 



(14) 
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The randomised approximation scheme is said to be a fully polynomial randomised 
approximation scheme, or FPRAS, if it runs in time bounded by a polynomial in 
\x\ and e~ l . Note that the quantity 3/4 in Equation (JTJJ could be changed to any 
value in the open interval (|, 1) without changing the set of problems that have 
randomised approximation schemes [221 Lemma 6.1]. 

We say that a real number z is efficiently approximate if there is an FPRAS for 
the problem which maps any input to the output z. Approximations to real num- 
bers are useful. For example, if q and 7 are approximations to q and 7 satisfying 



e n+mq < q < e n+mq 

and 



g n+m ry < < g n+m 

and 7 e = 7 for every e E E then 

(15) e ^Tutte ^Tutte 

Thus, to approximate ZTutte(G r ; g, 7), it suffices to first compute rational approxi- 
mations q and 7, and then approximate Z Tutte (G; g, 7). 

When the parameters are efficiently approximable reals, it is possible to approxi- 



mate quantities associated with the gadgets r and r that we studied in Section I4T31 
and EH 



Lemma 12. Suppose q > 2 is an efficiently computable real. Consider the gadget 
r from Section 4-4 with parameters t, N and g where g G [0, 1] is a rational 



number and N 1 ^ is an integer. There is an FPRAS for computing Z l and Z l , 
given inputs t, N and g. 

Proof. In this proof only, the notation rjv,t is used to make explicit the size of the 
gadget r. 

Let E(r Njt ) be the set of edges of -T/v,t- Let A k ' e denote the set of edge sub- 
sets A C E(r^ tt ) with k connected components containing vertices in T, and 
I other connected components. Let w(t,N,k,£) = Yl,A&A k < 1 Thus Z k = 

J2f =0 w(t,N,k,£)q e . We exhibit recurrence relations for w(t,N,j,k), from which 
it follows that these can be computed in polynomial time by dynamic program- 
ming. The rest of this proof is straightforward, but provides the details. To reduce 
the number of boundary cases we need to consider, it is convenient to allow one 
or other of k and t to take on the value —1. Of course, we stipulate 

w(t, N, -1, £) = w(t, N, k, -1) = 0, for all t, N,k,£> 0. 

Another easy-to-verify boundary case is 



w{t,0,k,£) 



when k = t and £ = 
otherwise, 
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which is valid for all t, k, £ > 0. Also w(t,N,k,£) = if exactly one of t and k 
is and w(0,N,0,0) is 1 is iV = and otherwise. Finally, w(0,0,k,£) = if 
k + £ > 0. 

The general recurrence, covering all situations other than the boundary cases 
already mentioned can now be given. 

(16) w(t,N,k,£) = ()( N ~i)w{i,j,l,0)w{t-i,N-j,k-l,£) 

l<i<t ^ ' 

l<j<N 

(17) + ( N ~l)w(0,j,0,l)w(t,N-j,k,£-l). 

(Although valid for all k,£ > 0, the recurrence becomes trivial when k + £ = 1, a 
point we must return to at the end.) The explanation is as follows. We partition the 
sum defining w(t, N, k, £) according to the connected component C in A containing 
some distinguished vertex v in K. (Note that iV > 0, since one of the boundary 
cases covers N = 0.) Let % = \C D T\ be the number of vertices of C that lie 
in T, and j — \C D K\ the number that lie in K. Summation ffT6]) deals with the 
situation i > and (fTT|) with the situation i = 0. The binomial coefficients in fTTol) 
count the number of connected components C 3 v that contain the distinguished 
vertex, and have the correct intersections \C D T\ = i and |C D if| = j with T 
and the factor w(i,j, 1,0) counts the weight of connected subgraphs of r[C}; 
and w(t — i, N — j, k — 1,£) the weight of subgraphs of r[K U T \ C] that have 
the correct number of connected components (i.e., k — 1 with vertices in T, and £ 
without). The analysis of summation (ITT)) is entirely similar. 

Assuming k + £ > 1, the recurrence is well founded, in the sense that all occur- 
rences of w(t', N', •, •) on the right hand side (apart from those that get multiplied 
by zero) have t' + N' < t + N. When k + £ = 1, i.e., (k,£) G {(0, 1), (1,0)} 
the recurrence (I16I17P becomes the trivial w(t, N, k, £) = w(t,N,k,£). Note that 
(k,£) = (0,1) entails t = (otherwise w{t,N,k,£) = 0) and (k,£) = (1,0) entails 
t > 0, so only one of the possibilities (k, £) G {(0, 1), (1, 0)} occurs for a given pair 
(t, N). In order to make progress in this situation, we apply a preprocessing step: 

w(t,N, 1,0)= (! + 7e)- w(t,N,k,i). 

e£E(F Nt ) 0<k<t,0<l<N 

k+l>l 

Exactly the same formal expression applies to w(t, N, 0, 1). The identity merely 
expresses complementation: the weight of subgraphs with k + £ = 1 equals the 
total weight of subgraphs, less the weight of subgraphs with k + t > 1. With this 
modification, the recurrence becomes well founded. 

Since only a polynomial number of distinct tuples (t, N, k, £) arise, the recur- 
rence can be solved by dynamic programming in polynomial time. Note that the 



20 



LESLIE ANN GOLDBERG AND MARK JERRUM 



computation of YleeE(r Nt )(^ + 7e) m ^ ne preprocessing step can be done exactly 
since g is a rational and iV 1//4 is an integer. 

Once the w(t, N, k, £) values are computed, we wish to estimate 

TV 

Z k = Y,w(t,N,k,£)q e . 
1=0 

Let e be the desired accuracy in the approximation-preserving reduction. Compute 
a rational q in the range e~ e ^ N q < q < e £ ^ N q and return Z k = ^2f =0 u)(t, N, k, £)q e . 

□ 

It will also be necessary for us to approximate the critical edge probability g, so 
that, with this approximation, the graph r approximately satisfies Equation (J7J) 
in Lemma O The following lemma shows that this is possible. 

Lemma 13. Suppose q > 2 is an efficiently computable real. Fix 7 > and let 
A = A c + (q — A c )/2. Suppose that x > is rational. Consider the gadget r from 
Section pJT5| with parameters t, N, and g. There is a randomised algorithm whose 
running time is at most a polynomial in x~ X > iV and t which takes input N and 
t (where it is assumed that N 1 ^ is an integer and that N > max{t 16 , iVo} for the 
constant N from Lemma\f^) and, with probability at least 3/4, computes a rational 
g in the range [iV~ 3 ,A/iV] such that, if A is drawn from RC(.T; q,p) , then 

(18) e- 7 < Pr ' y ^» - l -l < e * 7 
K ' 7 - Pi(Y(A) = i) ~ r 

Proof. First, we establish a useful preliminary fact. If < < 1/4, 0<<5<1 and 
i^sQ ^ Q ^ (1 + $)q then 

(19) e' 25 ^- < -t- < e 25 6 



1 — g 1 — g 1 — g 

To see this, note that e s g < 3/4 so 1 — e s g > g. Then letting x = e s — 1, 

l-g<l-g + x{l- e s g) - xg = (1 + x){l - (1 + x)g) = e\l - e s g), 

which, together with 1 + 5 < e s , gives the right-most inequality in fll9p . Similarly, 

ge~ s < 1/4 so 1 — e~ s g > g so letting x = 1 — e~ s , 

l-g>l-g + xg-x{l- e~ 5 g) = (1 - ar)(l - (1 - x)g) = e~\l - e~ s g), 

which, together with e~ s < j^, gives the left-most inequality in (fT9~|) . 

Now consider the gadget r with parameters t, N, and g. Let n and m be the 
number of vertices, and edges, respectively, of the variant r. We will be interested 
in three quantities which depend upon g: Z 1 , Z l and ( = ^y^w) • For this proof 
only, we will denote these as Z 1 , Z l Q and ( s to make the dependence on g clear. By 
Lemma CH C, Q = Z\jZ\. 
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Lemma H2] provides an FPRAS for computing ( g for rational parameters g. We 
will set 8 — x/ (16(n+m)) and we will use the FPRAS with accuracy parameter x/4 
to estimate ( g for every rational 

where fx is a non- negative integer with ™ (1 + 5) M < 4. Note that the number 
of integers in this range is 0(8^ \ogN). We will first "power up" the success 
probability of the FPRAS using standard techniques [21J so that the probability 
that any of these 0(8~ 1 log N) calls to the FPRAS fails is at most 3/4. If we find 
a rational g for which our estimate Q satisfies e _X//2 7 < ( g < e x ^ 2 7, we will return 
this value g. Clearly, this will be an acceptable answer, since e~ x ^( g < ( g < e x ^( g . 

To finish the proof, we just need to argue that we will try a g which we accept. 
Lemma M guarantees that there is a value g* G [iV~ 3 , XN^ 1 } such that ( g * = 7. In 
our computations, we will compute ( g for some g in the range j^g* < g < (1+6) q*. 
Recall that every edge e of the gadget has edge-weight 7 e = p(e)/(l — p(e)), and 
that p(e) = g for edges within the part of the gadget corresponding to the clique K . 
The weights of the other edges do not depend on g. Thus, by ffT§j) and by analogy 
to ffl"5l) . e~ x / 4 C, Q * < C, g < e x ^Cg*- Thus, g will be accepted. 

□ 



6. Approximation-preserving reductions and #BIS 

Our main tool for understanding the relative difficulty of approximation counting 
problems is approximation-preserving reductions. We use Dyer, Goldberg, Green- 
hill and Jerrum's notion of approximation-preserving reduction [S]. Suppose that 
/ and g are functions from S* to BL An "approximation-preserving reduction" 
from / to g gives a way to turn an FPRAS for g into an FPRAS for /. Here is 
the definition. An approximation-preserving reduction from / to g is a randomised 
algorithm A for computing / using an oracle for g. The algorithm A takes as input 
a pair (x,e) G S* x (0,1), and satisfies the following three conditions: (i) every 
oracle call made by A is of the form (w, 8), where w G X* is an instance of g, and 
< 8 < 1 is an error bound satisfying 8~ l < poly(|x|, (ii) the algorithm A 
meets the specification for being a randomised approximation scheme for / (as de- 
scribed above) whenever the oracle meets the specification for being a randomised 
approximation scheme for g; and (iii) the run-time of A is polynomial in \x\ and 
e-\ 

If an approximation-preserving reduction from / to g exists we write / <ap g, 
and say that / is AP-reducible to g. Note that if / <ap g and g has an FPRAS 
then / has an FPRAS. (The definition of AP-reduction was chosen to make this 
true). If / <ap g and g <ap / then we say that / and g are AP-interreducible, 
and write / =ap g- 
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The definitions allow us to construct approximation-preserving reductions be- 
tween problems / and g with real parameters without insisting that the param- 
eters themselves be efficiently approximable. Nevertheless, some of our results 
restrict attention to efficiently approximable parameters. According to the defini- 
tion, approximation-preserving reductions may use randomisation. Nevertheless, 
the reductions that we present in this paper are deterministic, except for where 
they make use of an FPRAS to approximate a real parameter. A word of warning 
about terminology: Subsequent to [8J, the notation <ap has been used to denote a 
different type of approximation-preserving reduction which applies to optimisation 
problems. We will not study optimisation problems in this paper, so hopefully this 
will not cause confusion. 

Dyer et al. [8] studied counting problems in #P and identified three classes of 
counting problems that are interreducible under approximation-preserving reduc- 
tions. The first class, containing the problems that admit an FPRAS, are trivially 
AP-interreducible since all the work can be embedded into the reduction (which 
declines to use the oracle). The second class is the set of problems that are AP- 
interreducible with ^Sat, the problem of counting satisfying assignments to a 
Boolean formula in CNF. Zuckerman [30J has shown that t^Sat cannot have an 
FPRAS unless RP = NP. The same is obviously true of any problem to which 
#Sat is AP-reducible. 

The third class appears to be of intermediate complexity. It contains all of 
the counting problems expressible in a certain logically-defined complexity class. 
Typical complete problems include counting the downsets in a partially ordered 
set [5J, computing the partition function of the ferromagnetic Ising model with 
varying interaction energies and local external magnetic fields [TT] and counting 
the independent sets in a bipartite graph, which is defined as follows. 

Problem: #BIS. 

Instance: A bipartite graph B. 

Output: The number of independent sets in B. 

We showed in [8] that #BIS is complete for the logically-defined complexity 
class #RHIli with respect to approximation-preserving reductions. We presume 
that there is no FPRAS for #BIS, but this is not known. 

7. Our results 

The next three sections of this paper give an approximation-preserving reduction 
from #BIS to the problem TUTTE(g, 7). 

We start by defining the problem of computing the Tutte polynomial of a uniform 
hypergraph H with fixed positive edge weights. For fixed positive real numbers q 
and 7 the problem UNlFORMHYPERTuTTE(g, 7) is defined as follows. 

Problem: UNIFORMHYPERTUTTE(g, 7). 
Instance: A uniform hypergraph H = (V, £). 
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Output: Z^ utte (H; q, 7), where 7 is the constant function with 7j = 7 for 
every / 6 E. 

In Section |H]we show that, for every q > 1, there is an approximation-preserving 
reduction from #BIS to UNIFORMHYPERTUTTE(g, q - 1). 

UNlFORMHYPERTuTTE(g, 7) is the problem of computing the Tutte polyno- 
mial of a uniform hypergraph. The most difficult part of the paper is reduc- 
ing this approximation problem to the problem of approximately computing the 
(multivariate) Tutte polynomial of an undirected graph (a 2-uniform hypergraph). 
Section |9] shows that for any positive real numbers q > 2 and 7 > 0, there is 
an approximation-preserving reduction from UNIFORMHYPERTUTTE(g, 7) to the 
following problem. 

Problem: Two WEIGHTFERROTUTTE(g) . 

Instance: Graph G = (V, E) with an edge-weight function 7' : E — > {7', 7"} 

where 7' and 7" are rationals in the interval [|V|~ 3 , 1]. 
Output: Z Tuttc (G; q, 7'). 

Finally, Section [10] gives an approximation-preserving reduction from the prob- 
lem TwoWEIGHTFERROTUTTE(g) to the problem TUTTE(g, 7) for any q > 2 
and 7 > 0. 

8. Approximately computing the Tutte polynomial of a uniform 

hypergraph 

We first define a parameterised version of t^BIS, and also a restricted version of 
this in which vertices on the right-hand side are required to have the same degree. 

Problem: #BIS(/x). 
Instance: Bipartite graph B. 

Output: Zis(B; /1) = J2i A 4 ' 7 ') where the sum is over all independent sets I 
in B. 

Problem: #SEMlREGULARBIS(/i). 

Instance: Bipartite graph B = (U,V, E) in which every vertex in V has the 
same degree. 

Output: Zis(B; fi) = A*' 7 '; where the sum is over all independent sets I 
in B. ^ 

Lemma 14. Suppose that > is efficiently approximable. Then #BIS <ap 

#SEMlREGULARBIS(/i) . 

Proof. In [HI Lemma 15], we gave an AP-reduction from #BIS to t^Bipar- 
titeMaxIS, the problem of counting maximum independent sets in a bipartite 
graph. Here, we give two AP-reductions - - first, a reduction from t^Bipar- 
titeMaxIS to 7^BIS(/i), and then a reduction from the intermediate problem 
#BIS(/i) to #SemiRegularBISGu). 
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Let B be an instance of ^BipartiteMaxIS with n vertices and m edges and 
let e be the desired accuracy of the approximation-preserving reduction. Let £ be 
the size of a maximum independent set of B and let Y be the number of maximum 
independent sets. To do the construction, we first need a rough estimate of /i, 
so compute a rational value /2 in the range |// < /2 < Let s be an integer 
satisfying 

n + 3 



1 < 



lg(l + 2/2/3) 



< s. 



Note that s > lg ^^ 2 ) • ^' ^ e ^ ne S ra P n w ith vertex set {(it, z) | u G 1^(5), z G 
[s]} and edge set {((«, i), (v,j)) \ (u,v) G E(B)}. Next, we need a more accurate 
estimate of [L. Compute a rational number fi in the range e - £ /( 60ns )^ < p, < 
e e/(60ns)^_ New since there are ns vertices in B', 

e-^°Z ls (B';^ < Z m (B';fi) < e^°Z ls (B'; /i). 

Thus, by using our oracle for #BIS(/i) with accuracy parameter e/60, we can 
compute a value Z satisfying 

(20) e-^ 30 Z ls (B'; jl) < Z < e^Z^B'; jl). 

Now note that every independent set of B' points out an independent set of B. 
The vertex u of B is in the independent set of B if there is at least one vertex (u, i) 
that is in the independent set of B' . A size-fc independent set of B thus makes 
a contribution of ((1 + [if — l) k to Zis(B'; ft). Since the number of independent 
sets of B is at most 2 n , we have 

(21) Z ls (B';fi) 2" <Y + - 

{ ] " ((i + A) s - i) c " ((i + A) s - i) " + 4' 

where the final inequality follows from the definition of s since fx > fi/2. The 
following simple procedure now gives a sufficiently accurate estimate to Y. Take 
the value Z from Equation (l2"Uj) . divide it by ((1 + jl) s — if and round down to the 
nearest integer. The fact that the accuracy is sufficient follows from (I2"0l and (|21[) . 
See Theorem 3] 

Finally, we present an AP-reduction from #BIS(/i) to ^SemiRegularBIS (/i) . 
Let B = (U,V,E) be an n-vertex instance of #BIS(/i) in which the maximum 
degree of a vertex in V is d > 1. We'll construct B' = (£/', V, E') - an instance 
of ^SemiRegularBIS(/x) in which every vertex in V has degree d. Let e be the 
desired accuracy of the approximation-preserving reduction. As above, compute a 
rational value fi such that |/i < \i < 1/2. Let /i~ denote the computed lower bound 

|/2 and let /x + denote the computed upper bound |/2. Let D{x) = (l + xf^ 1 , 
U(x) = xD(x), L(s,x) = (1 + x) s , and Y(s,x) = L(s,x) + D(x) — 1. We start by 
computing an integer s such that 

(22) 1 max( J D(/i+) - 1, U{^ + )) < 



L(s,pi~) 1 ' 6<in 
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Note that s = O^e" 1 ) so we can efficiently compute such a value s by starting 
with s — 1, and increasing s one-by-one until we find a value for which (|22p holds. 

Now let $ be a complete bipartite graph with vertex sets {z±, . . . , z^} and 
{y±, . . . ,y s }- For each vertex v G V of degree 5, take d — 5 new copies of 
and attach t> to vertex z\ of each copy. Let g < dn be the number of copies of \I/ 
that get included in 5'. 

Now Zis(\I / ; /i) = L(s, fi) + (l+fi)D(^) — l and the total contribution to Zis(^; fi) 
from independent sets including the vertex Z\ is U (/i). Let Y = Z IS (\I/; fi)-U(fi) = 
y(s,/i). Clearly, K is the total contribution to Z IS (\I/;/i) from independent sets 
not including the vertex Z\ 

Then 

Zis(B;fi)Y° < Z^B'-n) < Z ls (B; n)Z ls {% /i) 9 



so 



Zis(-B; a*) < — 77- — < ^is(-B; a*) 



Now 

• Zis(B'; /i) may be estimated with accuracy parameter e/3 using the oracle. 

• y 9 may be estimated directly (to the same accuracy) by computing a 
value ji satisfying e~ £ ^ 6dns ^n < ft < e £ /( Mns )/i which ensures that 

e -*/(6*0 L ( SjM ) < L( S) /t) < e £/(Mn) L(s,/i) 
And noting from fl22|) that 

L(s,ii) <Y = L(s,fi) (l + y ~M < L(s,//)e e/(6dn) . 



9 

< e £ / 3 , 



L(s,/i) 
Finally, ([22]) gives 

Z IS (*;//)V A ^(/i) 



which finishes the reduction. 



□ 



We are interested in the situation // > 1 so that we have g = /i+l>2in the 
parameters of UniformHyperTutte((?, 7) in Lemma fT5l below. 

Lemma 15. Suppose that fi > is efficiently approximate. Then 

#SEMlREGULARBIS(/!) <ap UNIFORMHYPERTUTTE(/i + 

Proof. Let B = (U,V,E) be an instance of #SemiRegularBIS(,u). Let H = 
(V, £ ) be an instance of UNlFORMHYPERTuTTE(/i + 1, /1) constructed as follows. 
Let s be a new vertex that is not in U U V and let V = U U {s}. For u 6 V, let 
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T(v) = {ueU : (u,v) G and let F v = T(v) U {s}. Let £ = LLy^- We wil1 
show below that 

(23) Z ls (B,fi) = (n + l) 

Now let e be the desired accuracy in the approximation-preserving reduction. 
To complete the reduction, we first compute a value x in the range 

e- £ / 2 Gu + l)- 1 < x < e e/2 (/i + I)" 1 . 

This is easy to do since [i is efficiently approximable and 

e s/2 l x e e/2 

< — 77: and 7- < 



fi + 1 ~~ e £ / 2 fi + 1 e~ £ / 2 fi + 1 _ jj, + 1 ' 

Then we use the oracle to estimate Zi utte (H; fx+1, /x) with accuracy parameter e/2. 

We finish the proof by establishing (|23|) . Let S C. V. The contribution to 
Zis(B,fj) from independent sets J with / fl V = S is ^l 5 l(// + l)l C7 l _ l r ( 5 )l ) where 
r(S) = U e5 I>)- That is, 



Z IS ( J B,/i) = ^^(/i + l 



ju|-|r(s)|_ 



On the other hand, the contribution to Z Tutte (H; /i + 1,/i) from the hyperedge 
set 8 = {F v : v G S} is 

^1(^ + 1)^ = ^(^ + 1)1^1-™+!^ 

since the vertices in r(S') together with s form one connected component, and all 
other vertices are isolated. Thus 



Z Tutte {H ]f i + l,fi) = ^/^(/i + l) 



v\-\r(S)\+i 
scv 

□ 



9. Approximately computing the multivariate Tutte polynomial of 

A GRAPH 

Lemma 16. Suppose that q > 2 and 7 > are efficiently approximable. Then 

UNIFORMHYPERTUTTE(g, 7) < A p TWOWEIGHTFERROTUTTE^) . 

Proof. Start with a t-uniform hypergraph H = (V, £). H is an instance of the 
problem UNlFORMHYPERTuTTE(g, 7). For convenience, let V = {vi, . . . ,v n } and 

£ = • • • ) /m}- 

The basic idea is to simulate each hyperedge with a copy of the gadget r from 
Section fl~3l (actually with the variant r from Section T4.4p . We start by setting 
the parameters. Let e be the desired accuracy in the approximation-preserving 
reduction (see the definition of an AP-reduction in Section [H]). Let \ = e/(4m) 
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The tolerance rj, which is a parameter to Lemma [TU1 may be chosen to be any 
value such that 

1 + ^ 1 + eX ^\ < e x, 

1-T] J 

Note that this inequality is achieved for rj = 0(e/m). 

The reduction will construct an instance G = (V, E) of the target problem 
TwoWElGHTFERRoTuTTE(g) along with an appropriate edge- weight function 

: E — > {7', 7"}. We will show that using an oracle to approximate the solution 
to this instance with accuracy parameter 6 = e/2 enables us to get within exp(±e) 
of Z Tuttc (#;g,7). 

Let A = A c + (q — A c )/2 and let N be the quantity from Lemmas M and [TUJ 
Let N be the smallest integer greater than max{t 16 , ^ _1//s , N } for which iV 1//4 is 
an integer. Using the algorithm from Lemma [T3] (suitably powered up so that 
its failure probability is at most 3/8) we can compute a rational g in the range 
[iV~ 3 ,A/iV] such that, if A is drawn from KC(r;q,p), then Equation (|18|) holds. 
Note from Lemma [10] that Equation dHJ) holds. 

The construction is as follows. For every j G [m], let Kj be a set of N vertices 
and let Vj = KjUfj. Let Ej = K^UKj x fa and let Gj be the graph (Vj, Ej). Note 
that Gj is a copy of the graph r from Section Wj$\ — fj is the set of terminals. Now 
we construct the graph G = (V, E) where V = V U UjsH anc ^ ^ = Ujg[ m ] 
Define 

if e G Kj 2 ^ for some j G [m], and 
if e g x for some j G [m]. 



p(e) 



Let 7 e = p(e)/(l — p(e)), and let 7' = {7 e } eg £;. Note that 7' = ^/(l — g) and 
7" = A^ _3 / 4 /(l — N~ 3 / 4 ). These are both rational values in the interval [| | 3 , 1], 
as required, since X/N < 1/2 by Lemma O 

Denote by IT 3 - the set of all partitions of fj. Given ttj G Hj, let A^ 3 be the 
collection of all edge subsets Aj C Ej that induce the partition 7Cj (into connected 
components) on fj. For A,- C let Y(Aj) denote the number of connected 
components that contain terminals (vertices in fj) in the graph (Vj,Aj). Let 
ft'(Vj, Aj) = n(Vj, Aj) — Y(Aj) be the number of remaining connected components 
(that do not contain terminals). For partitions it and tt', n V it' denotes the finest 
partition that is a common coarsening of 7r and tt'. If two elements are together 
in 7r or in tt' then they are together in the coursening. The "coarsest" partition 
has one block and the "finest" partition consists of singleton blocks. Technically, 
tt j G Hj is a partition of fj but we consider it as a partition of the entire vertex 
set V by extending it with singleton blocks. 

Now 

Ztuu c (G-<1^') = Y, ( 1 K{V,A) U^ 

ACE e£E 



K'(V h Aj) 
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Let Aj = A n Ej. Let j{Aj) = H eeA 7 e - Then 

^Tutte(G;g, 7 ')= E E g ^U, e[m] ^) JJ 7 (^; 

Vje[m] w " r J , 

where k{-k\ V • • • V 7r m ) denotes the number of blocks in the partition 7Ti V • • • V 7r m . 
Now let 

and pull out the contribution of each j to get 

m 

(24) *WG;g,Y)= E S^"^ H^)" 

TTj&lj j = l 

Vj6[m] 

At this point, we make some connections to Section WM Let Uj denote the set 
of all partitions of fj into k blocks and let Z* = XItt en fe Zji^j)- From Lemma [TT] 
we deduce that 

(o^ £l = F < Y ^ = k ) 

1 ; Z}' Pt(Y(A) = k') ' 

where A is drawn from RC(.T; q,p) for the function p from Section FOl 

Denote by _Lj the finest partition (with t blocks) of fj and by Tj the coarsest 
partition (with one block) on the same set. By Lemma [T3] Equation (lisp and 
Equation ( 1251) we have 



(26) e"*7 < ^4 < e* 7 . 

Use the FPRAS from Lemma [12] (again, suitably powered up so that its failure 
probability is at most 3/8) to compute a value a which satisfies 

(27) e~ x Zj(±j) <a< e*Z^). 

Now, define bj(0) = _L 3 - and 6 j ( 1 ) = Tj. Then, by restricting the sum in ([21)1 . 
to terms satisfying 7T,- G {JL 3 -, Tj} for all j, we get 

ZTutt e (G;q,~f')> E g K(fel(2l)v - vfcm(2m)) . 
ze{o,i} m 
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Letting ||z|| denote the Hamming weight of z and using Equation (|26p . the right- 
hand side is at least 

Z(_Lj) m e~ xm g K ^ 1 ^ 1 ^ v "' vbm(2m ^7" 2 ". 

z£{Q,l} m 



Using Equation (1271) . we get 

l(zi)V---Vft m (z m )) ||z| 

H 

z<E{0,l} m 



(28) Z Tuttc (G;q,~f') >a m e~ 2 * m ]T q<"^ - 



a 



m e- 2xm Z Tutte (H;q^). 



To finish, we just have to show that this lower-bound for ZT Uttc (G; q,j') is 
actually a good estimate because the terms that we threw away in inequality ( 1281) 
don't amount to much. For j G [m], let LTj = {T,} and let 11° = ILj \ {Tj}. 

Then, starting from 



in 



Z T utUG;qn , )= £ E g K(7riV - V7rm) n Z 



ze{o,i}™ eI # j=i 
Vj(z[m] 

Now note the partition ttj G II ■ j is refined by bj(zj). Thus, &i(zi) V ■ ■ ■ V b m (z n 
has at least as many connected components as tx\ V • ■ ■ V 7r m . So we get 

ZTutte(G;q,~f')< E £ ^^^-^^^n^^') 

Y?eH 



^ g K(6 1 ( Zl )v...v 6ro ( Zm )) ^ Y[Z j (7r J/ 



*e{o,i} m ^en^ i =1 

VjG[m] 
m 

^(fei(2i)v-v;. m ( 2m )) Zj(ii Zj ) 

zG{0,l} m j=l 



where Z,.(II*) = £ w . 6n ? z jM- Now, ^(Ilj) = ^(T,) < e 2 ^a. Also, ^(nj 
£* =2 z i • B y Equation ([25]), this is 



Pr(Y(A) = t) 
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where A is drawn from RC(i~';g,p) for the function p from Section 14.31 Since 
Zj < e x a, we have 

A >>- [ + Pr{Y{A) = t) 

Now by Lemma HO] Equation ©, J2t=2 Pv ( Y ( A ) = k ) < V- Thus > Pr(Y(A) = 
1) + Pt(Y(A) = t) > 1 - 7] and by Equation (0SJ, 

Pr(FM) = t) > 1 ~ V . 

Thus, 
Thus, 

ZTutte(G;g,7') < e 2 * m a m ]T g"(M*i)v-vM*»)) 7 IM| = e 2 * m a m Z Tutte (if ; g, 7 ). 

z€{0,l} m 

Thus, if we knew a quantity ?/> in the range 

e- £ / 2 < 1 < e £ / 2 , 

ZTuttelG; g,7') 

we would have 

~~ ^TutteCH"; 9, 7) ~~ 

so we have completed the AP-reduction. □ 

10. Approximately computing the Tutte polynomial of a graph 

In this section, we complete the approximation-preserving reduction from #BIS 
to TUTTE(g, 7) by giving an approximation-preserving reduction from the problem 
TwoWElGHTFERRoTuTTE(g) to the problem TuTTE(g, 7). 

For this, we first need to define what it means to "implement" an edge weight. 
This description is mainly taken from [131 Section 1.6]. Fix q > 2. Let W be a set 
of edge weights. For example, W might contain the edge weight 7 (a parameter 
of the problem TuTTE(g, 7)) or it might contain a selection of edge weights that 
we have already implemented using 7. Let 7* be a weight (which may not be 
in W) which we want to "implement". Suppose that there is a graph T, with 
distinguished vertices s and t and an edge-weight function 7 : E(T) — >■ W such 
that 

(29) 7* = qZst{T) 

where Z st (T) denotes the contribution to ZT U ttc(T; q, 7) arising from edge-sets A in 
which s and t are in the same component. That is, Z st (T) = ^ A •y(A)q K( - v ' A \ where 
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the sum is over subsets A C E(T) in which s and t are in the same component. 
Similarly, Z s \ t denotes the contribution to Z Tutte (Y; q,j) arising from edge-sets A 
in which s and t are in different components. In this case, we say that T and 7 
implement 7* (or even that W implements 7*). 

The purpose of "implementing" edge weights is this. Let G be a graph with 
edge- weight function 7. Let / be some edge of G with edge weight 7/ = 7*. 
Suppose that W implements 7*. Let T be a graph with distinguished vertices s 
and t with a weight function 7 satisfying (1291) . Construct the weighted graph G by 
replacing edge / with a copy of T (identify s with either endpoint of / (it doesn't 
matter which one) and identify t with the other endpoint of / and remove edge 
/). Define the weight function 7 as follows. 



So, as long as Z s \ t (T) is easy to evaluate and q is efficiently approximable, approxi- 
mating the multivariate Tutte polynomial of G with weight function 7 is essentially 
the same as approximating the multivariate Tutte polynomial of G with weight 
function 7. 

Two especially useful implementations are series and parallel compositions. 
These are explained in detail in [T7J Section 2.3]. So we will be brief here. Parallel 
composition is the case in which T consists of two parallel edges e\ and €2 with 
endpoints s and t and j ei = 71 and 7 e2 =72- It is easily checked from Equa- 
tion (j2UD that 7* = (1 + 7i)(l + 72) - 1. Also, the extra factor in Equation (13"0"|) 
cancels, so in this case Z Tuttc (G; q, 7) = Z Tuttc (G; q, 7). 

Series composition is the case in which T is a length-2 path from s to t con- 
sisting of edges e\ and t2 with j ei = 71 and 7 e2 =72. It is easily checked from 
Equation (j29~l) that w* = 7172/ (<7 + 7i + 72)- Also, the extra factor in Equation ( 130]) 
is q + 71 + 72, so in this case Z Tuttc (G; q, 7) = (q + 71 + 72)^Tutte(G'; q, 7). It is 
helpful to note that 7* satisfies 



Lemma 17. Suppose that q > 2 and 7 > are efficiently approximable. Then 

TWOWEIGHTFERROTUTTE(g) < A p TUTTE(g, 7). 

Proof. Let G = (V, E) be an instance of TwoWElGHTFERROTuTTE(g) with edge- 
weight function 7' : E — > {7', 7"} where 7' and 7" are rationals in the interval 




(30) ^T U t t e(G; q, 7) = —^Y-^-Z Tuttc (G; q, 7). 




We are now ready to prove this lemma. 
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[ | I 3 , 1]. We will assume without loss of generality that \E\ is sufficiently large 
with respect to the fixed parameters q and 7. 

let e be the desired accuracy in the approximation-preserving reduction. Let 
X = s/ (4(| Vj + \E\ 2 )) Let q be a rational in the range e~ x q < q < e x q and let 7 be 
a rational in the range e~ Xr j < 7 < e Xr y. Since q and 7 are efficiently approximable, 
the amount of time that it takes to compute q and 7 is at most a polynomial in 
\E\ and e" 1 . 

The idea of the proof is to show how to use series and parallel compositions 
from the set W = {7} to implement edge-weights 7'* and 7"* satisfying 

(31) e~ x i < 7'* < e x i 

and 

— y 11 ^ 11* ^ v " 

g *-ry < <y < g A 7 . 

Letting 7* be the edge-weight function derived from 7' by replacing 7' with 7'* 
and 7" with 7"*, Equation (fT5j) gives 

e-^TuttetG^y) < ^rutte(G; g,7*) < e £ / 4 Z Tutte (G; g, 7')- 
Let 7 be the edge-weight function which assigns every edge weight 7. We can 
think of our implementations as constructing a graph G such that Z^ ntte {G; q, 7*) 
is equal to the product of ZT U tte(G ! ; g, 7) and an easily-computed function of q 
and 7E We will ensure that each implementation uses at most \E\ edges, so the 
total number of edges in G is at most l-E"! 2 . To finish, we note (from (fT3|) ) that 

e" e/4 Z Tutte (L7; q, 7) < Z Tuttc (G; q, 7) < e e/4 Z Tutte (G; g, 7), 

where 7 is the constant edge-weight function which assigns every edge weight 7. 
We finish the approximation of Z^tteiG; q, 7) by using the oracle to approximate 
^Tutte(G; g, 7) using accuracy parameter 6 = e/2. 

It remains to show how to do the implementations. Taking 

We show how to use W = {7} to implement an edge-weight 7'* which satisfies 
7' — 7r < 7'* < 7'. This ensures that Equation (131 j) holds. The implementation of 
7"* is similar. 

Our implementation is taken from Section 2.1 of our paper [T3]. First, we can 
implement a weight 71 < 7 by taking a series composition of edges of weight 7 
for sufficiently large k. It suffices to take 

log(l + 4g) 



k 



log(l + g/7) 



1 This easily-computed function arises from the extra factor Z st (T) in Equation It is 

easy to compute because our implementations use only series and parallel composition 
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Then implement a weight 7j by taking a series composition of j copies of ji. The 
following (recursive) definitions are from [T5| Section 2.1] for integers j > 1. 



Then the implementation combines, in parallel, dj edges with edge- weight 7 j, 
for all j G [m] . 

The calculation in [T3"| Section 2.1] shows that the implemented value 7'* satisfies 
7' — 7r < 7'* < 7', as required. 

Now to finish we need to show that d\ H — • + d m < \E\, which we used above. 
First, note that the fixed parameters q and 7 give fixed upper and lower bounds 
on q and 7. Using the upper bound 7' < 1, we see that the value m is at most 
logarithmic in 7r _1 which is at most logarithmic in \V\, \E\, and e~ l . The calcula- 
tion in [T3"| Section 2.1] shows that the same is true of d\, . . . , d m . In fact, there is 
a fixed upper bound for dj depending only on q and 7. The proof of this fact uses 
the fact that < jj < 71 < ~. Since we assumed, without loss of generality, that 
\E\ is sufficiently large with respect to the fixed parameters q and 7, we conclude 
that di -\ — • + d m < \E\, as required. □ 

Proof of Theorem [IJ Theorem [T] follows from Lemmas HH [151 EH an d [13 D 



Lemmas and [13 have the following corollary. 

Corollary 18. Suppose that q > is efficiently approximable. Then #BIS <ap 
UniformHyperTutte^, q - 1). 

Thus, assuming that there is no FPRAS for #BIS, we can conclude that there 
is no FPRAS for computing the Tutte polynomial of a uniform hypergraph when 
the edge- weights are set to 7 = q — 1. The Ising model corresponds to the q = 2 
case of the Potts model. Thus, we conclude that there is no FPRAS for computing 
the partition function of the Ising model on a uniform hypergraph in which every 
edge has weight 1. 

We conclude this paper with a contrasting positive result for 3-uniform hyper- 
graphs. Consider the following problem. 

Problem: 3-UNIFORMHYPERTUTTE(g, 7) 
Instance: A 3-uniform hypergraph H = (V,S). 

Output: Z Tutte (iJ; q, 7), where 7 is the constant function with -fj- = 7 for 
every / G £. 

Lemma 19. Suppose that 7 > is efficiently approximable. There is an FPRAS 
for 3-UniformHyperTutte(2, 7) . 



iog((i + V)re:i 1 (i + 7.)^) 

log(l + 7i ) 



, and m 



log(g(l + V)/7r + l) 
log(g/7i + 1) 



11. 3-Uniform Hypergraphs 
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Proof. Jerrum and Sinclair [2U] have given an FPRAS for Tutte(2,7') for ev- 
ery 7' > 0. We will give a reduction from 3-UniformHyperTutte(2, 7) to 
Tutte(2, Y) where 7' = (1 + 7) 1/2 - 1. 

Let H = (V, £) be a 3-uniform hypergraph, an instance of the target problem 
3-UniformHyperTutte(2,7). Let 7 be the constant function with 7/ = 7 for 
every / G £. Let y — 7+ 1. Now, by Observation [21 

Z Tuttc (#; 2, 7) = Z Potts (H; 2, 7) = ^ ^monoM 

cr:V^{0,l} 

where mono (a) denotes the number of hyperedges / G 8 that are monochromatic 
in configuration a. 

Construct a (multi-)graph G with vertex set V and edge set 

E = IJ (V,w), (u,w)}. 

Let 7' be the constant function with 7^ = 7' for every j G -E. Let ?/' = 7' + 1 = 
y x l 2 . Now if a hyperedge / G £ is monochromatic in <r, it contributes y /3 to the 
corresponding term in Zp otts (G; 2, 7'). Otherwise, it contributes y' to the term. 
Thus, 

Z Tuttc (G; 2, 7) = Z Potts (L7; 2, 7') = y /|£:| Zp otts (#; 2, 7), 
which completes the proof. 

□ 
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